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Fractional (Levy-type) operators are known to be spatially nonlocal. This becomes an issue if 
confronted with a priori imposed exterior Dirichlet boundary data. We address spectral properties 
of the prototype example of the Cauchy operator (—A) 1 / 2 in the interval D = (—1,1) C R, with a 
focus on functional shapes of lowest eigenfunctions and their fall-off at the boundaries of D. New 
high accuracy formulas are deduced for approximate eigenfunctions. We analyze how their shape 
reproduction fidelity is correlated with the evaluation finesse of the corresponding eigenvalues. 


I. FRACTIONAL LAPLACIANS: R VERSUS D C R. 

The Fourier integral J R \k\^f (k)e~ lkx dk = —d^f(x)/d |x| M = |A|^/ 2 /(a;) is commonly interpreted as a definition 

of a fractional derivative of the /x-th order for p £ (0,2). The notation — (—A) M / 2 = — (A^/ 2 refers to a fractional 
Laplacian of order p/2 (restricted to dimension one i.e. to R) and two versions of a fractional dynamics (dimensional 
constants being scaled away): semigroup exp(—1|A| M / 2 ) / and unitary exp(— it\ A|^/ 2 ) /, 0-1- Here / stands for a 
Fourier transform of / € L 2 (R) and g{k) = \k\^f(k) is presumed to be L 2 (i?)-integrable. 

Apart from the unperturbed (free) case, the Fourier (multiplier) representation of the fractional dynamics has 
proved useful if an infinite or periodic support is admitted for functions in the domain, [2]. For the simplest quadratic 
(~ x 2 ) perturbation of the fractional Laplacian (the fractional oscillator problem), a complete analytic solution has 
been found in the specialized Cauchy oscillator case 3], by resorting to Fourier space methods. 

For more complicated perturbations, and likewise for a deceivingly simple problem of the fractional Laplacian in a 
bounded (spatial) domain, standard Fourier techniques seem to be of a doubtful or limited use, Q. A fully-fledged 
spatially nonlocal definition of the fractional Laplacian appears to be better suited to handle such problems, [§]-[l5|. 
See e.g. also [l6[ for a construction of Cauchy semigroups which arise from various perturbations of the Cauchy 
operator by bounded or locally bounded positive functions (i.e. external potentials). 


A. lAI' 1 / 2 on R. 


The fractional Laplacian — |Al^/ 2 , p £ (0,2) is a pseudo-differential (integral) operator and its action on a function 
from the L 2 {R) domain is defined as follows: 


-|A r / 2 /(x) = [ \f(x + y)-f(x)-^^}^(dy), 

Jr 1 + V 


( 1 ) 


where v^dx) stands for the Levy measure. This definition is commonly simplified by employing the Cauchy principal 
value of the involved integral (evaluated relative to the singular points of integrands) 


lAr'Vw = = 


dz 


Here, the Levy measure du^ has been made explicit and we point out a change of the integration variable y 
The Fourier representation of the the integral formula ((2]) takes the form 


iAr/M=- r < 1+ri8l, '“ r 


/ °° (p-iky _ 

^-i-TiX- 

-oo M 1+/i 


1 )dy 


The integral over dy, presuming its existence (which is not the case for p = 1) can be directly evaluated 


( 2 ) 


- X. 


( 3 ) 


( 4 ) 


Since Eq. © is undoubtedly valid for all p £ (0,1) U (1, 2), we can substitute back an outcome of © and employ an 
identity r(l+/z)r(— p) = — 7r/ sin(7r^t) (remember that the function T(—p) has simple poles at 0, —1, —2). Accordingly, 
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two potentially divergent entries compensate each other and the limit /r —»• 1 is now legitimate. Thus \A\^' 2 f(x) = 
jj= \k\^ f (k)e~ lkx dk for all fi £ (0, 2), and |fc| M is a Fourier multiplier of | A^/ 2 on R , as anticipated, see [HJ and 

The fractional Laplacian |A| M / 2 extends to a self-adjoint operator in L 2 (R) and induces a strongly continuous 
contraction semigroup exp(— 1\ A| M / 2 ) whose Fourier multiplier equals exp(—f |fc|^). 


B. Interlude: —A in D C R. 

The Hamiltonian-type expression H = —A + V, with V(x) = 0 for x £ D = (—1,1) C R, is an encoding of the 
Laplacian with the the Dirichlet boundary conditions (so-called zero exterior condition on R \ D) imposed on L 2 (R) 
functions f(x) in the domain of H: f(x) = 0 for \x\ > I. The problem is that so defined operator H , if restricted to 
a domain containing solely functions / £ L 2 (R) with a support in D, is not a self-adjoint operator in L 2 (R) (to this 
end we need to admit C^°(R) as a proper domain). 

We recall that C^°(R) comprises infinitely differentiable functions that are compactly supported in R. The notation 
Cq°(D) refers to a definite choice of the support to be D C R. The differential operator — A when acting in Cfi°(D) 
(we keep D = (—1,1)) defines a symmetric operator in L 2 (D). The problem of self-adjoint extensions in this case 
is a classic, c.f. [lf|. One of them, with a domain D(H) = {AC 2 [— 1,1],/(— 1) = 0 = /(1)} corresponds to a 
standard (quantum mechanical) infinite well problem; the AC 2 notation refers to an absolute continuity of / which 
gives meaning to the second derivative of / (of relevance at the boundary points of D). 

The spectral solution gives rise to the L 2 (D ) orthonormal eigenbasis, composed of real functions f n (x), n = 1,2,... 
such that f(x) = 0 for |x| > 1. More explicitly: f n (x) = cos(mrx/2) for n odd and sin(mrx/2) for n even, while 
respective eigenvalues read E n = (mr/2) 2 . It is clear that any f £ L 2 (D), in the domain of the infinite well 
Hamiltonian, may be represented as f(x) = c nfn(x), e.g. in the form of (trigonometric) Fourier series. 

At this point it is useful to mention that, quite independently of the self-adjointness issue, in L 2 (D) we have two 
inequivalent ways of making the Fourier analysis. If R\D is neglected and L 2 (D) is considered on its own (without any 
reference to L 2 (R)), then we can employ the previously mentioned Fourier series (e.g. the infinite well trigonometric 
eigenbasis). 

We shall pay more attention to the alternative approach. Namely, if L 2 (D) is considered as a subspace of L 2 (R), 
then for any 0 ^ f £ L 2 (D) we know that / £ L 2 (R) but / no longer belongs to L 2 (D). 

In fact, for any / £ Cq°(R) its Fourier transform / is an entire function (analytic on the whole complex plane). It 
does not vanish anywhere in R , except at infinity. If / would vanish on R\D, then necessarily it would be vanishing 
on D as well. 

Since / is infinitely differentiable it follows that / goes to zero along the real axis faster than the inverse of 
any polynomial. These properties hold true irrespective of the choice of the compact supporting interval D C R. 
Accordingly, TL 2 (D) C L 2 (R), but TL 2 (D) fj L 2 (D) = 0. 

Under the infinite well conditions, \k\ 2 still retains some of the Fourier multiplier (of —A) features. This view 
is supported by an approximation of the infinite well problem by a sequence of deepening finite wells, [l8i |. The 
convergence can be quantified in the L 2 (R) norm. 

Useful examples worked out in Bi, at the first glance, indicate an undoubtful relevance of k and k 2 multipliers. 
However, an the emergence of technical problems becomes conspicuous, if the multiplier property is to be elevated to 
eigenfunctions which are not Cq°(D). 

In fact, for the infinite well eigenfunctions f n (x) (L 2 (U)-normalized), their L 2 (R) Fourier transforms f n (k) = 
(2n)~ 1/2 f D f n (x) exp(— ikx) dx can be directly evaluated. We have ( f,f)L 2 (D ) = ( f,f)L 2 (R ) and there holds 
f R \fn(k )\ 2 k 2 dk = ( mr/2 ) 2 = (/„, (-A/„)). Likewise (/„, (-iV/ n )) = f R |/n(fc)| 2 kdk = 0 for all n> 1. 

One should be aware that the existence of mean values of Fourier multipliers does not mean that we can execute 
an inverse Fourier transform of e.g. the function g(k) = \k\ 2 f and retrieve —A f(x) as an image-function in L 2 (D). 
Actually, the inverse Fourier transform does not exist in this case, unless we adopt a weaker definition E3- We note 
that —A nf(x), at the boundaries ±1 of D, needs to be interpreted as a generalized function (distribution), and has 
a meaning only if smoothed out by a suitable test function. Incidentally, the eigenfunctions f n appear to play such a 
smoothing role, while evaluating mean values. 


C. Cauchy operator in D C R: from L 2 (R ) to L 2 (D). 

The previously indicated jeopardies related to the Fourier multiplier definition (in case of g = 2) surely extend to 
fractional Laplacians. Therefore, in the presence of spatial restrictions upon domains of nonlocal operators, we choose 
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to investigate their properties directly in configuration space with no recourse to Fourier transforms (and multipliers). 

The action of the Cauchy operator on C^ D (R) functions (differentiable, with continuous derivatives of all orders 
and compactly supported) is given by a specialized version of Eq. d2j): 


|A| 1 / 2 /(x) 


1 j- f(x) - f(x + z) d _ = 1 f f(x) - f(z) 


dz, x £ R. 


7r 


7T 


(5) 


and clearly has a Fourier representation with the multiplier |fc| f(k). We note in passing that so defined (—A) 1 / 2 
extends to an unbounded self-adjoint operator in L 2 (R). 

The Cauchy operator | A| x / 2 if restricted to a domain comprising solely L 2 (R) functions with a support in D and 
vanishing on R\D is not a self-adjoint operator in L 2 (R). However, if we consider the action of | A| x / 2 on test functions 
/ £ Cq°(D), then the restriction | A|p 2 / of | A I 1 / 2 to D is interpreted as the Cauchy operator with the zero (Dirichlet) 
exterior condition on R\D and is known to extend to a self-adjoint operator in L 2 (D), [§|. 

The passage from Cfi°(R) to Cfi°(D) ultimately amounts to disregarding any R\D contribution implicit in the 
formal definition © and makes the usage the Fourier multiplier representation either clumsy or redundant. 

Let us consider the D versus R\D interplay in more detail, by considering the action of |A| 1//2 on these Cfi°(R) 
functions which are actually supported in D. The major problem we wish to address is an explicit spatial form of the 
eigenvalue problem for |A|^ 2 , interpreted as |A|^ 2 f = E f where E £ R + is an eigenvalue and / £ L 2 (D). No closed 
analytic solutions are here available and various approximate methods were invented to optimize approximations of 
’’true” eigenvalues and shapes of respective ’’true” eigenfunctions, Q-O- 

Each known to date approximate eigenfunction xp(x), a inua, by construction belongs to the domain of |A|7 
and obeys exterior Dirichlet boundary data. However, generically |A|^ 2 'i/> £ L 2 (D ) no longer respects those data, 
while such a property is definitely expected from an acceptable L 2 {D) approximation of the right-hand -side of the 
eigenvalue formula |A|^ 2 / = E f. 

This problem is typically bypassed in the mathematical literature, where one considers the spectral problem of 
finding an eigenfunction (or its optimal approximation) in a weaker form. This (somewhat relaxed) approach to the 
spectral problem stems from the adopted definition of the action of |A|d in its Cfi°(D) domain. Namely, one demands 
that for / in the L 2 (P ) domain of |A|d there exists |A|.d / £ L 2 (D) such that ( g , |A|d f)L 2 (D) = (|A| g, f)L 2 (D) holds 
true for g £ Cg°(D), 0,111. 

Let us tentatively consider the action of |A| x / 2 on Cq°(R) functions xp(x), supported in D = (—1,1). Accordingly, 
for all x £ (—1,1) we have: 


2 xp(x) 
A d i/j(x) = -- 


1—X 


- xp(x + y) 


7T 1 — X 2 


dy. 


( 6 ) 


' —1—X 


The integral in © should be understood as the Cauchy principal value evaluated with respect to 0, according to 

/-;-* = linwo + ■ 

Let us change the integration variable y = t — x in Eq. ©. We have: 


A d i/j(x) = 


2 ip{x) 
7T 1 — x 2 


/-1 


j>{x) - j>{t) 
[t — x) 2 


dt 


(7) 


where the R \ D and D contributions are now clearly isolated. The Cauchy principal value of the integral in Eq. 
0 is no longer evaluated with respect to 0, but with respect to x. The integral expression in Eq. 0 which is now 
restricted to t £ (—1,1) while x £ (—1,1) by definition, is a special (Cauchy) case of so-called regional fractional 
Laplacian for D, 0-[ll|. 


II. TRIAL ANALYTIC FORMS OF APPROXIMATE EIGENFUNCTIONS. 

A. The ground state function. 

In the present paper, we shall employ Eq. © as the definition of the Cauchy operator in action on functions 
xp £ 'D{D) 1 in the spatially restricted domain D = (—1,1). We no longer require ip to belong to C^°(D) 7 we need 
however xp(x) to be infinitely differentiable in D and identically vanish for \x\ > 1. Let us address the eigenvalue 
problem Afjip(x) = Exp(x), whose approximate solution will be sought for subsequently. 
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The term ^ in Eq. © is a sum of the geometric series ^(1 + x 2 + x 4 + ...), (—1,1). Therefore it seems 

natural to assume that a solution </>(x) of Adi/j(x) = Etjj{x) might be represented in the form of the power series 

OO 

i/)(x) = c nX n as well, with a proviso that i/j(x) needs to be identically 0 at the boundaries ±1 of D. Thus we allow 

n— 0 

ip(x) to have a domain D = [—1,1], 

__ OO 

The ground state is an even concave function, [9|,[lCj, therefore we actually expect that ip( x ) = X) c 2 n x2n - Since 

n —0 

Y>(±1) = 0, there holds cq + C 2 + ... = 0. Inserting i/) to the eigenvalue formula and keeping x € (—1, 1) we get 


i T ^'(x)z-r( x )i 

n J z 2 


r'{x) Ir-... 


dz H— 


7T 1 — X' 


; 1 p(x) = Elp{x). 


( 8 ) 


presuming that the integral (Cauchy principal value) and series summation can be interchanged, next setting x = 0, 
we arrive at the series expansion which defines the ground state eigenvalue E, given ip(x): 



With cq 0, we have 


[l-l 

(t 

+ vy 

+ ? + ■• 


Co 

V l 

3 

5 

/ 


(9) 


( 10 ) 


The series converge, which follows (via the D’Alembert criterion) 
x e [-1,1]. 


OO 

from the convergence of c 2 nX 2n 

n —0 


ip(x) for all 




o 



FIG. 1: Approximate ground states for Cauchy wells. Numbers refer to: 1- infinite well tp(x) of Eg. (Hill : 2, 3, 4 - finite wells 
with depths Vo = 5000, 10000, 20000 respectively, [l5j; 5 - infinite well proposal of In the right panel, the curve 5 is out of 
the frame. 


By independent arguments, we know that the ground state function should be close (loosely speaking) to the cos(a;) 
while away from the boundaries of D and (1 — x 2 ) 1 / 2 in the close vicinity of the boundaries, [2, 8j, [l5|. Let us consider 
the trial approximation of ground state function, given in an analytic form: 

i/>(x) = C\f(l — x 2 ) cos(aa:), (11) 

where the coefficient a has been adjusted to differ slightly from § — f-, known to be the leading term in the asymptotic 
eigenvalue formula for the infinite Cauchy well, 0: 

1443 /7T7r\7r n tt tt tt 

a = 4096^ V2 _ 8/64~256~512~ 1024 ~ 4096’ 


( 12 ) 
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C = 0.921749 being the L 2 (D) normalization constant. 

In the recent paper El we have introduced an algorithm for evaluating approximate eigenfunctions of finite Cauchy 
wells of arbitrary depth. The idea was to implement as close approximation of the infinite well spectral properties 
in terms of those for very deep finite wells. In Fig. 1 we have depicted approximate ground state functions for finite 
wells of depths Vq = 5000,10000, 20000, an approximate ground state for the infinite Cauchy well proposed in Ref. 
[8j and our approximate formula i/j(x) of Eq. (fill) for the infinite well. 

In the left panel, curves for Vo = 5000,10000, 20000 and that for if>(x) are graphically indistinguishable in the 
adopted scales, while the proposal of Q is conspicuously different. In the right panel a vicinity of the maximum has 
been enlarged and the proposal of Q, in the adopted scale, is out of frame. 

OO 

Let us expand the approximate ground state if of Eq. CD into power series ip(x) = ]T] C 2 n x 2n with x £ [—1,1]. 

n —0 

The expansion coefficients can be explicitly identified and we reproduce numerical values for first few of them: 

c 0 = 0.921749, c 2 =-0.743145, c 4 = 0.011510, c 6 =-0.020710, c 8 = -0.015567, . , 

c 10 = -0.012318, ci 2 = 0.009969, c 14 = -0.008234, ci 6 = -0.006922, ci 8 = -0.005910. 

Although if(x) is not a ’’true” eigenfunction but an approximation of the ground state, by employing merely first 10 
expansion coefficients in the series expansion Eq. m, we obtain the very rough outcome E = 1.15318. The ground 
state eigenvalues have been approximated by other (more accurate) methods and, up to four decimal digits we have 
e.g.: 1.1577 according to 0, 1.1573 according to E3- 

We point out that a convergence of the series CD is very slow. To get more accurate approximation of the eigenvalue 
associated with the approximate ground state ^}( x ) (fill) we need to account for much larger number of expansion 
coefficients c 2 n. 


B. Analysis of Anip{ x )- 


At the moment we are not that much interested in producing high accuracy approximate formulas (this issue will 
be addressed subsequently in the present paper). The analytic expression Eq. (flD for ip(x) is extremely useful for 
another purpose. 

Namely, we can make explicit the action od Ad upon functions with definite geometric shapes and analyze not only 
how much Ajjif(x) deviates from i^{x) and ultimately from Eip(x) (with a properly adjusted value E ), but also the 
boundary behavior of those functions. See e.g. [12i Il3| for some hints in this connection. 

Given tp(x) = Cy/ (1 — x 2 ) cos(ax), |x| < 1, we would like to know whether the Dirichlet boundary data (e.g. 
vanishing of a function for \x\ > 1) are respected by Ad ip(x)- To this end let us consider 

— X+l 

1 f 1 p(x) — 1 p(x + z) 

7T J Z 2 

—x—1 —x—1 


l f if{x) - cy 1 - (x + z) 2 {l - 7 z(x + z) 2 - y 4 (x + z) A - .. .) 
7 T J Z 2 


where we have expanded y/cosa(x + z) into power series whose coefficients are denoted 72 ra , (70 = 1). 
First few coefficients are given explicitly, 


a 2 a 4 19a 6 559a 8 

72 - Y’ 74 ~ %’ 76 ~ 5760’ 78 ~ 645120' ^ ^ 

We integrate each term of the series separately. The integral corresponding to 70 = 1 can be rewritten as follows (see 
e.g. HU) 


— X + l _ 

ib(x) f 1 — \/p + qz + rz 2 

- / -- 2 - dz 

7i ■ J z 2 

—x—1 


(16) 


where p = 1/cos(aa:), q = —2x/{l — x 2 ) cos(ax), r = —1/(1 — x 2 ) cos(a:c). We evaluate the integral in the sense of 
its Cauchy principal value (see e.g. [23|), temporarily skipping the factor ip(x)/Tr: 


— X + l 


—x—1 


1 ~ Vp + qz + rz 2 _ _ 2 


(1 — x 2 ) y/(l — x 2 ) cos(aa;) 


(17) 
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We note that the first term in the above (after restoring ip/ 71 ) cancels its negative in the defining expression ([ 6 ]) for 
A D tp{x). 

Subsequent integrals can be evaluated analogously, but with ip{x) fully incorporated in the integrated expressions. 
We merely disregard (but keep in mind) an omnipresent coefficient C/n and denote a = 1 — x 2 ,b = — 2 x,c = —1. 
With this proviso other integrals follow: 


72 


-x+l __ 1 

/' \Ja + bz + cz 2 (x + z) 2 f u 2 

/ -79-= 72 


\/l — u 


(it — x) : 


2 7r(l — 6 a; 2 ) 

du = 72 ---. 


(18) 


—ai+1 


I Va + bz + cz 2 (x + z) 4 r " ,4 

74 / -79- = 74 


~ U 2 


du = ~- 7 r(l + 12a ; 2 — 40a; 4 ), 


(it — a ;) 2 


(19) 


— tc + 1 


, Va + bz + cz 2 (x + z) 6 r "A 

76 / - 5 -= 76 


U”V 1 — U 2 , 76 


. du = -^ 7 r(l + 6 a ; 2 + 40a; 4 — 112a; 6 ). 

(it — a ;) 2 16 

— X— 1 —1 

Accordingly, after reintroducing the factor C/tt we arrive at the polynomial expansion of A^ip 

OO 

A D 1 p(x) = C , y^ 72 nW 2 n(x), 

n=0 


( 20 ) 


( 21 ) 


where coefficients 72 n for 0 < n < 4 have been explicitly identified before while W 2 n {x) are polynomials of degree 2n, 
like e.g. 

. 1 — 6a; 2 . . 1 + 12a; 2 — 40a; 4 . 1 + 6x 2 + 40a; 4 — 112x 6 

wo{x) = 1, w 2 {x) = - - -, w 4 ( x) =---, w 6 {x) = -—-• (22) 





°-To ^05 OO 05 bo 
X 



FIG. 2: Left panel: a comparison of i/j(x ) = CV (1 ~ a; 2 )cos(aa;) (dotted line) and An'i’ (solid line). Middle panel: \Ad^{x) — 
Eip(x) | with E = 1.156. Right panel: supremum of | Ad ip — E tp(x)\(a) for E = 1.156. The a-axis is scaled in units 7 t/4096. 


The series expansion of \/ cos(aa;) converges very fast for x £ [—1,1]. Accordingly, by accounting for only first few ex¬ 
pansion terms we get quite good approximation of Ad V- I n Fig. 2 , we have depicted both ip{ x ) = Cy/(1 — x 2 )cos(ax) 
(a dotted line) and the resultant Ad V (solid line). 

Since in Ref. ITfl we have numerically identified an approximate ground state eigenvalue to be close to E = 1.156 
(eventually correctable to E = 1.1573, 0 ), let us employ the latter value instead of less accurate E = 1.15318 
obtained in a rough reasoning presented before. Then, we can point-wise compare V(x) against AdV( x ) by depicting 
a curve | AdV( x ) ~ Eijj(x)\,x G D, see e.g. Fig. 3. 
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The deviation of A^ipi^x) form Eip(x) appears to be rather small and effectively concentrates in the vicinity of the 
boundaries x = ±1. Since we have lim A d i/j( x) = 0.130753, lim Eif>(x) = 0, there holds \A d i/j(x) — Eip(x)\ < 

x —>-±1 : r —>-±1 

0.130753, x £ D which is the best point-wise estimate ever obtained in the literature on the (shape) subject, compare 
e.g. 0 (Lemma 1, formulas 8.9 and 8.10) and [gj. We emphasize that the behavior of \Adi/j(x) — Eif >| is fairly 
robust with respect to the specific choice of E. The dominant contribution comes to the upper bound comes form the 
behavior of Ad if){x) at the boundaries od D. 

The shape of the approximate ground state, as proposed in [ 10 , while away from the boundaries (i.e. around 
x = 0) is significantly different from our present finding and from the numerically-deduced behavior of eigenfunctions 
in finite but deep Cauchy wells, [Tsj]. 




FIG. 3: Left panel: i/j(x) is represented by dotted lines (black and red variants are practically indistinguishable), Ao'tpix) for 
a = 14437t/ 4096 (solid black) and a = 150l7r/4096 (solid red). Right panel: \Ad^>{x) — Et/)(x)\ for E = 1.156 and previous as, 
respectively in black and red. 


We note that in the definition (11111 we can in principle vary a. Taking a supremum of \Ad ^(x) — Eip(x)\(a) over 
x £ D as a criterion for how close A^ip is to Eip, we realize that the optimal a choice would be a = 150l7r/4096, see 
e.g. at the location of the minimum in Fig. 4. 

To see better how the choice of a may affect the shape of Api^ix) and \Ad^{x) — Eip(x) |, we display the behavior 
of these functions comparatively for a = 14437t/ 4096 and a = 150l7r/4096. We note that for a = 150l7r/4096 
we have \Ad'4>{x) — Eijj(x)\ < 0.06 which is much better point-wise estimate than previously obtained 0.13 (for 
a = 14437t/ 4096). The price paid is slightly worse fitting away form the ±1 boundaries. 


C. First excited state. 

The ground state function, previously denoted ip(x), in fact should be labeled by n = 1, hence denoted il>i{x). To 
avoid notational confusion, the first excited state (n = 2) will be denoted ^(x). 

We introduce a trial analytic approximation of a ’’true” first excited state in the form: 

ip 2 {x) = — C sin(/3x) \J{\ — x 2 ) cos(/3x), (23) 


where 


1760 7T 7T 7T 

^ = 4096 7 ’’ = 2 _ 16 ~ 128’ ^ ^ 

and C = 1.99693 is a normalization constant. The minus sign is basically irrelevant, but is introduced for graphical 
comparison purposes. The parameter [3 fitting comes from our data for deep but finite Cauchy wells, [l5j . In Fig. 
6 we display various approximate formulas for the ’’true” excited eigenfunction, comparing our analytic guess with 
approximate curves proposed in Q and [lB}. 













FIG. 4: Approximate expressions for the first excited state function. Numbers refer to: I - ip 2 (x) (1231) . 2 - finite Cauchy well 
outcome, depth Vo = 5000, [l5|], 3 - an approximation proposed in Q. Left panel: magnified vicinity of a maxium. Right panel: 
magnified vicinity of the boundary x = 1 . 


D. Analysis of Ad ip2(x). 

Main arguments follow these of Section II.A and II.B. The first excited function ip 2 (x) is odd, hence its power 

OO 

series expansion i/) 2 (x) = c 2 n+ix 2n+1 contains only odd-labeled coefficients c n . Like in section IIA, we deduce the 

n —0 

eigenvalue E 2 


1 - — 

(t 


C 7 

+ -r + . . 

)} 

Cl 

\ 1 

3 

5 

) 


E 2 = 


Since we know explicitly the numerical values of coefficients C 2 , let us list few of them 


ci = —2.695662, 
cn = 0.010042, 


c 3 = 3.394555, 
ci 3 = 0.008479, 


c 5 = -1.040718, 
C 15 = 0.007567, 


c 7 = 0.152499, 
ci 7 = 0.006774, 


c 9 = 0.000755, 
cig = 0.006095, 


(25) 


(26) 


and next insert them directly to the expansion (1261) . while disregarding the remainder of the series. An approximate 
eigenvalue reads E 2 = 2.72874, to be compared (even though the result is very rough) with the deep finite Cauchy 
well prediction E 2 = 2.7534 of Ref. [H} (eventually correctable by 0.0013 to E 2 = 2.7547) and E 2 = 2.7547 of Ref. 
Q. The series (l26l) converge slowly, therefore much larger number of coefficients c n need to be accounted for, to make 
reliable the numerical value for i? 2 . 

To deduce Adi[> 2 (x), let us first analyze the integral 


— fC+1 


1p2 {x) -4> 2 (X + Z) 


dz 


— X—l 


alone, see e.g. © for comparison. For clarity, in the power series expansion 

OO 

sin(/3x)v / cos(/3a;) = ^ 72n+ia; 2ri+1 , x e D, 

n=0 

we enlist few 72 n+i in their explicit numerical form: 

l 3 19/3 5 

77 = 


. 5 /? 3 19/3 5 

7i ~ ^ 73 ~ 12 ’ 75 _ 480 ’ 


31/3 7 


8064 


(27) 


(28) 


(29) 


The Cauchy principal value can be evaluated for each expansion term of sin(/3a;)ycos {fix) separately. Accordingly, 
for the first term we have 


—ai+l\ 


— C 71 lim 
£—>-0 


cVl — x 2 — (x + z) 1/1 — (x + z ) 2 


dz. 


\—x—l 


(30) 
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Since 


— £ —CC+1\ 


cVl — • 


-dz = xyjl - x 2 ( —-- ^ + 


1 — x 2 


(31) 


and 


— £ — X+l\ 


~(x + z)y/l ~ (x + z) 2 = ( y/1 - ( X-e ) 2 + y/i - (x + £) 2 


V— x —1 


+ 2nx + u(x, e), (32) 


where (here unspecified) u(x, e) approaches 0 if e —► 0, for all x € D. 

The first term in (1511) . if multiplied by —C^i/ir, cancels its negative in the expansion of 2ip2(x)/n(l — x 2 ). Because 
of 


lim 

£—>-0 


2-\/l — x 2 — yjl — (x — e) 2 — y/1 — (x + e) 2 


= 0, 


(33) 


the first expansion term has an ultimate form —Cjiwi(x), where w\{x) = 2x. 

In connection with (1531) we point out that troublesome (divergent) 2/e entries (related to the Cauchy principal value 
evaluation) are cancelled by their negatives coming from the principal value procedure of the form J5]) while adopted 

~ X+1 (x+z) ~ X+1 ib ( ) 

to — f ' f ’ 2 /T z dz. The remaining expansion terms of f V2 )f' dz are cancelled by their negatives that originate 

— X —1 — X—1 

from i/j 2 (x)/(l — x 2 ). 





FIG. 5: Left panel: a comparison of the approximate function i/> 2 (x) (dotted line) with (solid line). Middle panel: we 

depict \Ad ^2 — E 2 ip 2 {x)\, where E = 2.75. Right apnel: supremum of \Ad^I> 2 {x) — E 2 tp 2 {x)\, E = 2.75 as a function of /?. The 
horizontal axis is scaled in units 7r/4096. 


An ultimate form of A Dll’ 2 {x) is 


Ad^x) = ^2 l2n+lW2n+li.x), X £ D, 


(34) 


n —0 


where the coefficients 72 n+i, n = 0,1,... have been introduced before (while expanding sin(/?a;)-\/cos(/3a:)) and first 
few polynomials w 2 n+i have the form: 


wi(x) = 2x, W 3 (x) = —x(l — 4x 2 ), w^{x) = 


c(— 1 — 8x 2 + 24a: 4 ) 


, w 7 (x) = 


;(—1 — 4a; 2 — 24a; 4 + 64a; 6 ) 


■ (35) 


The convergence of (1341) is much worse than that encountered in connection with the ground state function. There¬ 
fore a number of polynomials employed in the approximation of ip 2 must be relatively large to make that approximation 
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reliable. In Fig. 7, we compare Ajj^ix) with ip 2 (x), for an approximation restricted to first 15 series expansion terms 
only. 

Like in case of the ground state function, we ask for an affinity of Ajyip 2 (x) with F^V^x), where we adopt the value 
E 2 = 2.75. In Fig. 8 the affinity function \Ad4>2(x) — E 2 ip 2 (x]\ is depicted and found to be bounded point-wise by 
0.1462 which is much better estimate than any ever obtained, ma¬ 
in the definition of an approximate function ip 2 ( x ) we have still some flexibility allowed with respect to the choice of 
the parameter /?. In Fig. 9, sup \Ae>4’2(x) — E 2 ip 2 {x)\ is depicted as a function of /3 in the vicinity of /3 = 17607r/4096. 
A minimum is achieved for /3 = 1762-7T/4096 and sets the upper bound 0.1344. 

In principle, we can provide analytic approximations for (consecutive) higher excited state functions. However, our 
discussion of Section II should be considered merely as a warm-up, a preparatory step to address more serious goals. 

Let us note that all ultimate formulas have involved the polynomial expansions. Interestingly, we could not associate 
them with any members of a hypergeometric family or other orthogonal polynomials in D = (— 1,1)). We shall follow 
this polynomial expansion strategy in the next section to get most accurate to date approximations of eigenvalues 
and eigenfunctions in the infinite Cauchy well problem. 

We shall impose one more demand, deliberately absent in the existing literature on approximate Cauchy well 
eigenfunctions. We need that actually not only ip{ x ) is to vanish identically for |x| > 1, but Adi/j(x) as well, to 
become in all respects as close as possible to E'i/j(x), where E stands for an approximate eigenvalue. 


III. POLYNOMIAL EXPANSIONS OF EIGENFUNCTIONS IN THE INFINITE CAUCHY WELL: 

PUSHING AHEAD APPROXIMATION FINESSE. 

A. Ground state function addressed anew. 

We wish to solve the eigenvalue problem Adi^(x) = Eip(x), with the Cauchy operator Ad in D 1 defined in Section 
II. This means that in the serch for approximate solutions, an approximation accuracy can be made arbitrarily fine, 
with the growth of the degree of the truncated polynomial expansion of the sought for eigenfunction. Functions in 
the domain of Ad are restricted by the exterior Dirichlet condition ip{x) = 0, |x| ^ I, but we impose the very same 
restriction upon the resultant Ad t/>(x), given ?/;(x) £ D. 

We take Eq. © as a working definition of Ad and proceed with its integral part, here denoted 


— X + l 

B D ip{x) = i J - 


l/)(x) — ljj(x + z) 


dz. 


(36) 


— X— 1 

Let us consider the action of Bd upon functions ip{ x ) = x 2n \/l — x 2 . We get: 

r 2 


BdV 1-x 2 = 


+ 


BdX 2 \J\ — x 2 = — 


7r 1 — x 2 

2 x 2 \J\ — x 2 1 — 6x 2 
7T 1 — x 2 2 


i-, 4 A- o 2 x 4 \/l — x 2 1 + 12x 2 — 40x 4 

Bdx V1 — x z = - 


BdX^\J 1 — x 2 = — 


7T 1 X 2 8 

2 xVl - x 2 1 + 6x 2 + 40x 4 - 112x 6 


7T 1 — X 2 16 

Accordingly we have 

/ - 2 \/\ 

B D x 2n Vl - X 2 = - ; — - T -1- (c 2 „ + 3c 2 n-2X 2 + 5c 2rl -4 + ... + (2n + 1 )c 0 x 2n ), 


7r 1 — x 2 

where C 2 n are expansion coefficients of the Taylor series for \/l — x 2 . Namely, we have 

" _(2 n)\ 

n - 

n=0 n=0 


(37) 

(38) 

(39) 

(40) 

(41) 


v^l - X 2 = C 2n x 2n = 


J2n 


(1 — 2n)(n!) 2 4 n ' 


(42) 
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which allows to rewrite Be)X 2n \/l — x 2 as follows 


Be)X 2n \/l — x 2 = 


2 x 2n Vl^x 2 ^ (2fc)!(2n + l-2fc) 2n _ 2k 
7r 1 — x 2 (1 — 2k){k\) 2 A k 


Then ground state function is even, hence we can expect its power series expansion in the form: 


(43) 


i/j(x) = C \/1 - x 2 Y ct 2 nX 2n , ao = l. (44) 

n —0 

where our major task is to deduce the expansion coefficients c^n- 

Coming back to the definition (6]), of Ae >, we realize that Ae^{x) = Be>^{x) + 2^(x)/7r(l — x 2 ). The second term 
of this expression, upon employing the trial i/j(x), as in Eq. (14T1) or (42), clearly cancels the first term of Be>'^{x), 
compare e.g. (36)-(42). 

In view of this, the action of Ae> upon the ground state candidate-function i/j(x) of (TH1) greatly simplifies. Ultimately, 
the eigenvalue problem Ae>'4’{x) = Eip(x) takes the form: 


v - ' (2fc)!(2n + 1 - 2k) 2n - 2 k 

to 2n h (l-2*)(*!) a 4* 




(2fc)! 


2k 


(1 — 2fc)(fc!) 2 4 /c ‘ 


E 

n —0 


^2 


which can be re-ordered as follows 


EE a2 " 

k— 0 n—k 


(2k)\(2n + 1 — 2k) 2n _ 2k 
(1 — 2fc)(fc!) 2 4 fc 


oo oo 


EE Ea. 2 n 

k —0 n —0 


{%k)\ 2k+2n 

(1 — 2fc)(fc!) 2 4 fe 


(45) 


(46) 


This system of equations (from which the coefficients a. 2 n are to follow for all n ) looks hopelessly discouraging, if we 
are interested in a fully-fledged solution of the eigenvalue problem. However, things simplify if we look for approximate 
solutions, readily accessible upon a truncation of otherwise infinite series. 

By definition we know that any solution ^(x) is defined in the domain D = [—1,1] and obeys the boundary condition 
z/)(± 1) = 0. We extend this restriction to Ae>i/j(x) and demand 

lim Ae>i/j{x) = 0. (47) 


In what follows we shall assume that ao = 1- This choice is possible, in view of the presumed normalization (we have 
C involved). Each considered truncated series approximating ip^x) will be l 2 (D) normalized, this operation being 
encoded in a multiplicative constant C. 


Example 1: Let us exemplify our procedure by truncating the series after the polynomial of degree 2. We have 
Ad{ 1 + ot 2 X 2 )\/l — x 2 = 1 — a .2 (^ — 3cc 2 ) and the condition (1-TT1) implies 1 + = 0. Accordingly we end up with an 

approximate eigenstate ip(x) = C( 1 — 2/5x 2 )v / l — x 2 where C = i/875/996 « 0.937291. The approximate eigenvalue 
reads E = 1.2. 

Example 2: An analogous procedure for series terminating at the polynomial of the 4-th degree gives rise to 
Ae>( 1 + a 2 x 2 + a 4 X 4 )Vl — x 2 and (1T71) implies 1 + |«2 + yr a 4 = 0. Moreover, we have 1 — \a .2 — ^ 0:4 = E and 
3«2 — |«4 = E [— ^ + 0 : 2 )- The coefficients readily follow with values 0:2 ~ —0.353189 and 014 ~ —0.0346746. 
The approximate eigenvalue reads E ^ 1.18093. The normalized approximate eigenfunction takes the form 
ip{x) = C( 1 - 0.353189x 2 - 0.0346746x 4 )Vl - x 2 , where C = 0.931331. 

It is clear, that the procedure can be continued indefinitely by increasing the polynomial degree at the series 
truncation ’’point”. An approximation accuracy grows with the polynomial degree. The polynomial degree growth 
increases the number of linear equations to solve. 

Let (ak, n ) denote a matrix with elements 


We recall that 


(277. T 1 2/c )Ck 


(2k)\(2n + 1 - 2k) 
(1 - 2k){k\) 2 A k ’ 


n ^ k. 


(2*Q! 

(1 - 2fc)(fc!) 2 4 fc ' 


(48) 


( 49 ) 
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If we consider an approximation of ip(x) by series terminating at the polynomial of degree 2n, the eigenvalue problem 
we address takes the form of a linear system of equations with unknown E and a 2 n (we recall our choice of a 0 = 1) : 

n i 

^ ( Oi2kO'k—i,k — E ^ ( OL2kCi—k 5 ^ — 0, 1, . . . , Tl 1, 

k=i fc=0 

n / m \ 

^ ( f 0^2m ^ I 0* (b0) 

m—0 \ fe=0 / 

The last identity in (15U1) comes from our demand (1T71) . here adopted to 1 — £ 2 = 0 at ie = ±1. 


- 

C 

E 

OL2 

OL4 

«6 

«8 

aio 

OU2 

a 14 

OL 16 

W2 

0.937291 

1.200000 

-0.400000 

- 

- 

- 

- 

- 

- 

- 

W4 

0.931331 

1.180929 

-0.353189 

-0.03467461 

- 

- 

- 

- 

- 

- 

Wq 

0.927253 

1.170127 

-0.333863 

-0.00891937 

-0.0332900 

- 

- 

- 

- 

- 

w s 

0.925363 

1.165443 

-0.326159 

-0.00332500 

-0.0173088 

-0.0221718 

- 

- 

- 

- 

mo 

0.924339 

1.162981 

-0.322268 

-0.00097523 

-0.0134732 

-0.0111668 

-0.0163303 

- 

- 

- 

U>12 

0.923728 

1.161534 

-0.320035 

0.00025555 

-0.0117497 

-0.0084661 

-0.0081667 

-0.0126748 

- 

- 

W14 

0.923337 

1.160614 

-0.318637 

0.00098488 

-0.0107978 

-0.0072137 

-0.0061348 

-0.0063114 

-0.0102120 

- 

wi 6 

0.923071 

1.159993 

-0.317704 

0.00145367 

-0.0102098 

-0.0065016 

-0.0051721 

-0.0047139 

-0.0050726 

-0.0084590 

Wl8 

0.922884 

1.159555 

-0.317051 

0.00177313 

-0.0098192 

-0.0060507 

-0.0046131 

-0.0039456 

-0.0037753 

-0.0041958 

W20 

0.922746 

1.159234 

-0.316576 

0.00200068 

-0.0095458 

-0.0057448 

-0.0042525 

-0.0034927 

-0.0031448 

-0.0031160 

W30 

0.922409 

1.158447 

-0.315422 

0.00253637 

-0.0089180 

-0.0050710 

-0.0035043 

-0.0026342 

-0.0021154 

-0.0017921 

W40 

0.922868 

1.158159 

-0.315006 

0.00272257 

-0.0087053 

-0.0048519 

-0.0032737 

-0.0023882 

-0.0018495 

-0.0015006 

W50 

0.922230 

1.158022 

-0.314810 

0.00280842 

-0.0086084 

-0.0047537 

-0.0031724 

-0.0022828 

-0.0017390 

-0.0013839 

W60 

0.922198 

1.157948 

-0.314703 

0.00285494 

-0.0085562 

-0.0047014 

-0.0031190 

-0.0022279 

-0.0016822 

-0.0013250 

W70 

0.922179 

1.157902 

-0.314638 

0.00288292 

-0.0085249 

-0.0046702 

-0.0030874 

-0.0021957 

-0.0016492 

-0.0012910 

w 80 

0.922166 

1.157872 

-0.314595 

0.00290106 

-0.0085047 

-0.0046501 

-0.0030671 

-0.0021751 

-0.0016283 

-0.0012696 

Woo 

0.922158 

1.157852 

-0.314566 

0.00291348 

-0.0084909 

-0.0046364 

-0.0030534 

-0.0021612 

-0.0016142 

-0.0012552 

WlOO 

0.922152 

1.157837 

-0.314545 

0.00292235 

-0.0084810 

-0.0046267 

-0.0030437 

-0.0021514 

-0.0016042 

-0.0012451 

W150 

0.922137 

1.157802 

-0.314496 

0.00294331 

-0.0084578 

-0.0046039 

-0.0030208 

-0.0021285 

-0.0015811 

-0.0012218 

W200 

0.922132 

1.157789 

-0.314478 

0.00295063 

-0.0084497 

-0.0045959 

-0.0030130 

-0.0021206 

-0.0015732 

-0.0012139 

W300 

0.922129 

1.157781 

-0.314466 

0.00295585 

-0.0084440 

-0.0045903 

-0.0030074 

-0.0021151 

-0.0015677 

-0.0012083 

U>400 

0.922127 

1.157778 

-0.314461 

0.00295767 

-0.0084419 

-0.0045884 

-0.0030055 

-0.0021132 

-0.0015658 

-0.0012064 

W500 

0.922127 

1.157776 

-0.314459 

0.00295851 

-0.0084410 

-0.0045875 

-0.0030046 

-0.0021123 

-0.0015649 

-0.0012056 


TABLE I: Approximate solutions of the eigenvalue equation Adi!>{x) = Ei/j(x). The approximating polynomial of degree 
2 n is indicated by W 2 n{x). We report first few values of coefficients ct 2 k for each 2n-th case, together with an approximate 
eigenvalue E and the normalization coefficient C. For comparison we report the ground state eigenvalue reported in Ref. [ill ]. 
E= 1.157773883697 (based on a diagonalization of the 900 x 900-matrix). Our ultimate result actually is E — 1.1577764. 


The system (l50l) can be solved for various series truncation choices, up to the 2 n = 500 polynomial degree. All 
computations have been carried out by employing the routines of Wolfram Mathematica, which appear to be dedicated 
to solving even very large linear systems of equations. 

Remark: One needs to be aware that m, as a system of n + 1 equations, has more than one solution. We select 
an optimal approximation of the ground state function by selecting a solution with then least value of E. The same 
system of equations produces solutions that approximate higher (excited) even eigenfunctions. There appear also 
complex solutions which we discard as physically irrelevant. 

Our findings are gathered in Table I, where we report explicit values for first few expansion coefficients ct 2 n of 
approximating polynomials, the approximate eigenvalue E and related normalization constant C. Symbols W 2 n refer 
to approximating polynomials of degree 2 n. The computed eigenvalues definitely drop down with the growth of the 
polynomial degree 2n, with a visible stabilization tendency. In our opinion, our result is much sharper than that 
















n 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

02n 

-0.3144595 

0.00295851 

-0.0084410 

-0.0045875 

-0.0030046 

-0.0021123 

-0.0015649 

-0.0012056 

-0.0009571 

-0.0007784 

020+2 n 

-0.0006454 

-0.0005439 

-0.0004647 

-0.0004016 

-0.0003506 

-0.0003088 

-0.0002740 

-0.0002449 

-0.0002201 

-0.0001990 

040+271 

-0.0001808 

-0.0001650 

-0.0001512 

-0.0001391 

-0.0001284 

-0.0001189 

-0.0001104 

-0.0001029 

-0.0000960 

-0.0000899 

O60+2n 

-0.0000843 

-0.0000793 

-0.0000747 

-0.0000705 

-0.0000666 

-0.0000631 

-0.0000598 

-0.0000568 

-0.0000540 

-0.0000515 

080+271 

-0.0000491 

-0.0000469 

-0.0000448 

-0.0000429 

-0.0000411 

-0.0000394 

-0.0000378 

-0.0000363 

-0.0000349 

-0.0000336 

O100+277. 

-0.0000324 

-0.0000312 

-0.0000301 

-0.0000291 

-0.0000281 

-0.0000272 

-0.0000263 

-0.0000255 

-0.0000247 

-0.0000239 

0 120+277 

-0.0000232 

-0.0000225 

-0.0000219 

-0.0000213 

-0.0000207 

-0.0000201 

-0.0000196 

-0.0000191 

-0.0000186 

-0.0000181 

0 140+277 

-0.0000177 

-0.0000172 

-0.0000168 

-0.0000164 

-0.0000160 

-0.0000157 

-0.0000153 

-0.0000150 

-0.0000147 

-0.0000143 

O160+2ti 

-0.0000140 

-0.0000137 

-0.0000135 

-0.0000132 

-0.0000129 

-0.0000127 

-0.0000125 

-0.0000122 

-0.0000120 

-0.0000118 

0 180+277 

-0.0000116 

-0.0000114 

-0.0000112 

-0.0000110 

-0.0000108 

-0.0000106 

-0.0000104 

-0.0000103 

-0.0000101 

-9.97* 10' 6 

0200+277 

-9.81 * 10 -6 

-9.67* 10' 6 

-9.53* 10“ 6 

-9.39* 10' 6 

-9.26* 10' 6 

-9.13* 10 -6 

-9.00* 10' 6 

-8.88* 10' 6 

-8.76* 10 -6 

-8.65* 10' 6 

0220+277 

—8.54 * 10 -6 

-8.43 * 10~ 6 

-8.33* 10” 6 

-8.23* 10' 6 

-8.13* 10' 6 

-8.03 * 10 -6 

-7.94 * 10~ 6 

-7.85 * 10~ 6 

-7.76* 10” 6 

-7.68 * 10~ 6 

0240+277 

-7.59* 10 -6 

-7.51 * 10 -6 

-7.43 * 10" 6 

-7.36 * 10" 6 

-7.28* 10~ 6 

-7.21 * 10" 6 

-7.14* 10 -6 

-7.08* 10~ 6 

-7.01 * 10" 6 

-6.95* 10" 6 

0260+277 

-6.88* 10' 6 

-6.82 * 10~ 6 

-6.77* 10 -6 

-6.71 * 10~ 6 

-6.65 * 10' 6 

-6.60* 10“ 6 

-6.55* 10' 6 

-6.50* 10' 6 

-6.45* 10 -6 

-6.40 * 10~ 6 

0280+277 

-6.35* 10' 6 

-6.31 * 10~ 6 

-6.27* 10 -6 

-6.22* 10' 6 

-6.18* 10' 6 

-6.14* 10 -6 

-6.10* 10' 6 

-6.07* 10~ 6 

-6.03* 10” 6 

-6.00 * 10~ 6 

0300+277 

-5.96* 10 -6 

-5.93 * 10" 6 

-5.90* 10" 6 

-5.87* 10~ 6 

-5.84* 10~ 6 

-5.81 * 10" 6 

-5.79 * 10~ 6 

-5.76* 10~ 6 

-5.73* 10" 6 

-5.71 * 10“ 6 

O320+277 

-5.69* 10' 6 

-5.67* 10' 6 

-5.65* 10“ 6 

-5.63 * 10~ 6 

-5.61 * 10~ 6 

-5.59* 10 -6 

-5.57* 10' 6 

-5.56* 10' 6 

-5.54* 10 -6 

-5.53* 10' 6 

0340+277 

-5.51 * 10~ 6 

-5.50* 10~ 6 

-5.49* 10 -6 

-5.48* 10' 6 

-5.47* 10' 6 

-5.46* 10 -6 

-5.45 * 10~ 6 

-5.45 * 10~ 6 

-5.44* 10 -6 

-5.44* 10' 6 

0360+277 

-5.43* 10 -6 

-5.43 * 10 -6 

-5.43* 10" 6 

-5.43* 10“ 6 

-5.43 * 10~ 6 

-5.43* 10" 6 

-5.43* 10' 6 

-5.44 * 10~ 6 

-5.44* 10" 6 

-5.45* 10 -6 

0380+277 

-5.46 * 10 -6 

-5.46* 10' 6 

-5.47* 10 -6 

-5.48 * 10~ 6 

-5.49 * 10~ 6 

-5.51 * 10 -6 

-5.52 * 10~ 6 

-5.54* 10' 6 

-5.56* 10 -6 

-5.57* 10~ 6 

0400+277 

—5.59 * 10~ 6 

-5.62 * 10~ 6 

-5.64* 10 -6 

-5.66* 10' 6 

-5.69* 10' 6 

-5.72* 10 -6 

-5.75 * 10~ 6 

-5.78* 10~ 6 

-5.82* 10 -6 

-5.86* 10' 6 

0420+277 

-5.90* 10 -6 

-5.94 * 10 -6 

-5.98* 10" 6 

-6.03* 10~ 6 

-6.08* 10' 6 

-6.13* 10" 6 

-6.19* 10 -6 

-6.25 * 10~ 6 

-6.32* 10" 6 

-6.39* 10" 6 

0440+277 

-6.46* 10 -6 

-6.54* 10' 6 

-6.62* 10 -6 

-6.71 * nr 6 

-6.81 * 10~ 6 

-6.91 * 10 -6 

-7.03* 10' 6 

-7.15 * 10~ 6 

-7.28* 10 -6 

-7.42 * 10~ 6 

0460+277 

-7.57* 10 -6 

-7.74* 10' 6 

-7.92* 10" 6 

-8.12 * 10~ 6 

-8.35 * 10~ 6 

-8.59* 10" 6 

-8.87* 10~ 6 

-9.18* 10~ 6 

-9.54* 10" 6 

-9.95 * 10~ 6 

0480+277 

-0.0000104 

-0.0000110 

-0.0000117 

-0.0000125 

-0.0000136 

-0.0000150 

-0.0000171 

-0.0000204 

-0.0000271 

-0.0000540 


TABLE II: For the approximate ground state function, the corresponding polynomial W500 is displayed in detail, in terms of its expansion coefficients 02 n - Note a 
numbering convention: in the first row we have displayed consecutively 02,04... up to 020- 


OJ 
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FIG. 6: Left panel: a comparative display of polynomials W 2 n(x) of degrees 2,4, 6, 10, 20, 30, 50, 70,100, 150, 200, 500 and the 
curve \J cos (14437ra;/4096) (gold) which has been a building block in the formula (1111) . Middle panel provides an enlargement 
in the vicinity of the right boundary. Right panel depicts various approximations of the ground state function at the right 
boundary x = 1: 1 - curve Cw 500 (*)%/! — x 2 , 2 - curve of Q, 3 - ipi{x) ~ (1 — l*!) 1 ^ 2 of 0, 4 - Vo = 500 finite well ground 
state of |15l| . 


reported in Ref. a see also Numerical values of the coefficients ct^n grow as well with 2 n growth, with a clear 
stabilization tendency. 

In Table II we make explicit the functional form of the polynomial rusoo- All coefficients 0 : 2 n are reproduced as well. 
At the moment that provides the best available approximation of the ground state function in the infinite Cauchy 
well problem. 

Since we are interested in fine details of the eigenfunction shape, it is instructive to display the behavior of the 
major eigenfunction building blocks, i.e. the polynomials W 2 n {x), in the vicinity of the boundaries of D. In Fig. 
6, we admit x from the exterior of D , i.e. |x| > 1. We compare the near-the-boundaries behavior of polynomials 
of degrees 2,4,6,10, 20, 30, 50, 70,100,150, 200, 500, with a function \J cos (14437nr/4096) (colored gold) appearing in 
the definition of the trial ground state function (fill) (c.f. Section II). That clearly explains deficiencies of the cosine 
factor in the adopted definition and obvious virtues of the present polynomial expansion method. We note that the 
polynomial degree growth, is accompanied by a steeper decent of the representative curves at the boundaries. 

At this point it is also instructive to have a comparison of the boundary behavior of the approximate ground state 
function proposed in the literature so far. We note that ip{x) of Ref. Q, at the boundaries, is fapp (for all practical 
purposes) identical with our ^>(x) ~ Cw 5 oo(x)y/l — x 2 . 

As before, we can analyze a deviation of Ajji/j(x) from Eip(x) in terms of | Ar>ip(x) — Eip(x )|, x £ D. Results are 
depicted in Fig. 7. Note that for the polynomial of degree 2 n = 500 we get an upper bound \Ajj 4’( x ) — Eip(x )| < 0.01. 
We note that with the growth of the polynomial degree, the near-the-boundary maximum of \Adi/j(x) — Eip(x)\ 
decreases. 


B. Other even eigenfunctions ip 2 k+\{x), k > 1. 

The system (50) of equations has been dedicated to obtain even eigenfunctions. As mentioned before it has infinitely 
many solutions, both real and complex-valued. Each real solution is interpreted as an approximation of a certain 
eigenfunction. Since for each resolved polynomial of degree 2n, we can jointly compute approximate eigenfunctions 
and the corresponding eigenvalues, there appears a natural ordering with respect to increasing E-values which we 
enumerate by consecutive odd numbers 2k, with k = 1 corresponding to the ground state. That allows for a systematic 
selection of higher rank even eigenfunctions. We keep intact the notation W 2 n(x) for an approximating polynomial of 
degree 2n, although one should keep in mind that for each consecutive E 2 k+i, we deal we the corresponding 2 k+ 1-th 
polynomial (and appropriate 2k + 1-th set of expansion coefficients). The previous Table II data refer to the ground 
state solution ipi(x) only. The coefficients tables for other polynomials area available upon request. 
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FIG. 7: Left panel: | Anip(x) — Eip{x )| where ip = Cy/1 — x 2 W 2 n(x), with 2 n = 2, 4, 6,10, 20, 30, 50, 70,100, 150, 200, 500. Right 
panel: polynomial degrees 2 n = 50, 70,100,150, 200, 500, ip(x) in the vicinity of the right boundary x = 1. 





FIG. 8: An approximation of the third eigenfunction ipz(x). Numbers refer to: 1 - ip 3 ~ Cwsoo{x)y/l — x 2 , 2 - ips(x) according 
to 3,3-finite Vo = 500 Cauchy well ground state, [TB]. Left panel: enlargement of the vicinity of the minimum. Right panel 
- enlargement of the vicinity of x = 1. 


C. Odd eigenfunctions, ip 2 k{x ), k > 1. 


In the notation of Section III.A, the odd eigenfunctions case i can be handled by invoking: 


D /I- 2 2 X ^ 1 ~ x2 , O 

Boxyl — ar =-;-^-1- 2x , 


7r 1 — x 2 


n S P, -9 2 X 3 Vl — X 2 ( 1 n , 

B D X 3 \J 1 - x 2 = --- 5 -b 2a; -- + 2a; 2 , 

7r 1 — ar V 2 


B D x 5 \/l - a; 2 = - - 7 y 1 + 2a; - 2 • ^a; 2 + 3a: 4 

7T 1 — ar V 8 2 




-l-2.^ 2 -3.ia; 4 + 4 x f 


(51) 

(52) 

(53) 


7r 1 — x 2 


H- 2x 


(54) 
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FIG. 9: \Ad'4 , 3(x) — Ez%j)z{x)\ where ip 3 is a product of C\J\ — x 2 and a polynomial of degree 2n, we depict 2 n = 
10, 20, 30, 50, 70,100,150, 200, 500. Right panel refers to 2 n = 70,100, 150, 200, 500. Note that for 2 n = 500, we have 
\A D tp3{x) — E3i/j 3 (x)\ < 0.07. 





FIG. 10: For an approximate eigenfunction ip 5 (x) we display: 1 - 2n = 500, 2 - ips^x) of Ref. @1, 3 - the fifth finite Vo = 500 
Cauchy well eigenfucntion (computed, but not reproduced in EH)- Left panel - minimum vicinity enlargement. Right panel - 
x = 1 vicinity enlargement. 


i.e. 


- 2 x 2n+1 \/l - x 2 

B D x 2n+1 \j 1 - X 2 =- ; --- 5 -b 2x(c 2n + 2c2n-2X 2 + 3 c 2n -4 + ... + («+ 1 )c 0 x 2n ), (55) 


7 r 1 — x 2 


where c 2n are Taylor series coefficients for y/1 - x 2 . We are interested in odd eigenfunctions and seek them in the 
form: 


'ipjx) =C\/l-x 2 y^/3 2 n+iz 2ra+1 , /3i = 1. 


(56) 


n —0 


As in the case of even functions, we look for solutions of the eigenvalue problem Ajjip{x) — Eip(x), so arriving at 


/Wi Y 

n—0 k =0 


( 2 fc)!( 2 n + 2 — 2k) 2ra _ 2fc+1 
(1 — 2fc)(fc!) 2 4 fe 


=e± 7 


„ 0 (i - 2km w „ 0 


Y fc n + ix 


2n+l 


(57) 
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FIG. 11: \Ad'4 ) 5{x) — E s ip s (x)\ where i/js^x) is inferred for 2n = 20,30,50,70,100,150,200,500. Right panel -2 n = 
70,100,150, 200, 500. 





FIG. 12: ip 2 (x). numbers refer to: 1 - polynomial of degree 501, 2 - according to @J, 3 - finite Cauchy well Vo = 500, 0- Left 
panel - enlargement of the minimum. Right panel - enlargement of the x = 1 boundary. 


and next at 


E E ^2«+ 1 

k—0 n—k 


(2fc)!(2?r + 2 2k) 2 n-2k+i 

(1 — 2fc)(fc!) 2 4 A ' 


EE E fcn+l 


(2 k)\ 


k—0 n —0 


(1 - 2/c)(fc!) 2 4 fc * 


2k+2n+l 


Additionally, we impose our boundary condition 


(58) 


lim AdiJj(x) = 0. (59) 

X —»±1 

The procedure adopted to find polynomial approximations of eigenfunctions and eigenvalues in the even case, can 
be extended to the odd case as well. Let 

„ „ IN (2fc)!(2n + 2 — 2k) , , N 

h, n - (2n + 2 - 2 k)c k - ^ _ 2fc)(A!) 2 4 fe ’ n ^ ( 60 ) 

where 

(2fc)! 

k (1-2fc)(fc!) 2 4 fe ' 


( 61 ) 
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FIG. 13: \Ad' 4 , 2 (x) — E 2 if> 2 (x)\, for polynomial degrees 2n + 1 = 11, 21, 31, 51, 71,101,151, 201, 501. Right panel - 2n + 1 = 
71,101,151,201,501. 



FIG. 14: 3 / 14 ( 3 :), numbers refer to: 1 - polynomial of degree 501, 2 - tp 4 (x) according to ||, 3 - finite Cauchy well Vo = 500 
outcome, 15]. Left panel - minimum enlargememnt. Right panel - right boundary enlargement. 


The linear system of equations from which all /? 2 n+i and E are to be inferred has the form (we set /3i = 1) 


n i 

^ ^ /^2fc+l bk—i,k = E ^ ^ (^2k-\-l^i—k 1 Z = 0, 1, . . . , 77- 1, 

k=i k =0 

n / m \ 

=0. (62) 

m=0 \ fc=0 / 

Like in case of (50), the present system of equations has infinitely many solutions, real and complex-valued for each 
fixed n. As before, in the set of real solutions, an ordering relation is set by referring to an increasing order of 
computed eigenvalues E 2 k k > 1. 

Skipping unnecessary repetitions of previous arguments, we display our findings concerning ip 2 {%) and to¬ 

gether with estimates for \Ajjip{x) — Eip{x)\. A comparison is made with previously reported results on the shape 
of corresponding approximate eigenfunctions. It appears that our method provides most accurate to date data 
for both approximate eigenfunctions and eigenvalues and provides the sharpest to date point-wise estimates for 
\A D ip(x) - Eip(x )|. 

For completeness, in Table III we report first five computed eigenvalues ordered against the approximating polyno- 
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°-T0 o!o 05 -L0 

X 


FIG. 15: \Au'4 ) a{x) — E 4 ip 4 (x)\ where polynomial degrees 
2n + 1 = 71,101,151,201,501. 



odd, 2n+ 1 = 21,31,51,71,101,151,201,501. Right panel: 


n 

Si 

e 2 

Ez 

e 4 

S 5 

7 

1.160614 

2.768252 

4.351150 

5.946117 

7.337136 

8 

1.159993 

2.765561 

4.344362 

5.934918 

7.584192 

9 

1.159555 

2.763594 

4.339381 

5.928041 

7.512343 

10 

1.159234 

2.762114 

4.335613 

5.922546 

7.509991 

15 

1.158447 

2.758299 

4.325845 

5.907535 

7.485347 

20 

1.158159 

2.756826 

4.322066 

5.901342 

7.475242 

25 

1.158022 

2.756110 

4.320233 

5.898233 

7.470144 

30 

1.157948 

2.755709 

4.319211 

5.896463 

7.467238 

35 

1.157902 

2.755463 

4.318584 

5.895363 

7.465432 

40 

1.157872 

2.755301 

4.318173 

5.894634 

7.464235 

45 

1.157852 

2.755188 

4.317889 

5.894127 

7.463403 

50 

1.157837 

2.755107 

4.317685 

5.893759 

7.462802 

75 

1.157802 

2.754913 

4.317196 

5.892875 

7.461356 

100 

1.157789 

2.754844 

4.317024 

5.892559 

7.460842 

150 

1.157781 

2.754795 

4.316900 

5.892331 

7.460473 

200 

1.157778 

2.754777 

4.316857 

5.892251 

7.460343 

250 

1.157776 

2.754769 

4.316837 

5.892214 

7.460282 

K 

1.157773 

2.754754 

4.316801 

5.892147 

7.460175 


TABLE III: We display the computed eigenvalues Ek, k > 1 under an assumption that polynomials of degree 2 n were 
employed in the definition of even eigenfunctions and 2n + 1 for odd eigenfunctions. The capital K in the last line indicates 
data collected from Ref. 0 - 


mial degree and compare them with those obtained by other arguments in Refs. S0. It is clear that by increasing 
the polynomial degree n we can achieve still higher finesse level of a computational accuracy with which both eigen¬ 
functions and eigenvalues can be retrieved. It is seen that with the growth n, there are definite stabilization symptoms 
in the numerical outcomes for the eigenvalues. We would like to point out that to increase the reproduction accuracy 
of higher ’’true” eigenfunctions, one should pass to higher than n = 500 polynomial degrees. The same pertains to 
the 0.01 upper bound for AdiJj(x) — Eip{x)\ if ip( x ) is an approximation of the ground state. To push that bound 
closer to 0, higher polynomiola degrees are necessary. 
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IV. OUTLOOK 

In the present paper we were largely motivated by: (i) on the one hand -successes in the approximate evaluation of 
eigenfunctions and eigenvalues for the infinite Cauchy well problem §|, fljl , (ii) various drawbacks in the physics- 
motivated procedures to solve that eigenvalue problem, summarized in [lj, |l5| . It has been often mentioned that the 
” true” eigenfunctions show a striking similarity to trigonometric sine or cosine functions (identified as eigenfunctions 
of the standard Laplacian in the interval) when away from the boundaries of D, while their fall-off towards zero at the 
boundaries should be similar to \/l — x 2 . Our trial function considerations of Section II proved that the trigonometric 
connection is somewhat deceiving, since the square root of a trigonometric function has been involved. In Section III 
we have resolved the away-from-the-boundary behavior by means of truncated polynomial expansions, that bear no 
obvious similarity, neither to trigonometric functions nor to their square roots. Things became more complicated and 
subtle, since the polynomial shapes actually dictate minute details of the eigenfunctions fall-off to 0 at the boundaries. 
In this connection, we point out the peculiar fall-off of approximating polynomials around ±1, as depicted in Fig. 6. 
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